Mon. Not. R. Astron. Soc. OOP. [TlfSTI C2002') Printed 20 January 2013 (MN I*TeX style file v2.2) 



The Tumultuous Formation of the Hubble Sequence at 
z > 1 Examined with HST/WFC3 Observations of the 
Hubble Ultra Deep Field 

C. J. Conselice^*, A.F.L. Bluck^'^, S. Ravindranath^, A. Mortlock^, A. Koekemoer 
F. Buitrago^, R. Griitzbauch^, S. J. Penny^'^ 

^ University of Nottingham, School of Physics & Astronomy, Nottingham, NG7 2RD UK 
^ Gemini Observatory, Hilo Hawaii 96720, USA 

^ lUCAA, Pune University Campus, Pune 4.11007, Maharashtra, India 
^Space Telescope Science Institute, Baltimore, MD USA 21218 

^Endeavour Fellow, Swinburne University of Technology, Hawthorn, Victoria 3122, Australia 
Accepted ; Received ; in original form 

ABSTRACT 

We examine in this paper a stellar mass selected sample of galaxies at 1 < z < 3 
within the Hubble Ultra Deep Field, utilising WFC3 imaging to study the rest-frame 
optical morphological distribution of galaxies at this epoch. We measure how appar- 
ent morphologies (disk, elliptical, peculiar) correlate with physical properties, such as 
quantitative structure and spectral- types. One primary result is that apparent mor- 
phology does not correlate strongly with stellar populations, nor with galaxy structure 
at this epoch, suggesting a chaotic formation history for Hubble types at z > 1. By 
using a locally defined definition of disk and elliptical galaxies based on structure and 
spectral- type, we find no true ellipticals at z > 2, and a fraction of 3.2 ± 2.3% at 
1.5 < 2 < 2. Local counterparts of disk galaxies are at a similar level of 7 — 10%, much 
lower than the 75% fraction at lower redshifts. We further compare WFC3 images with 
the rest-frame UV view of galaxies from ACS imaging, showing that galaxies imaged 
with ACS that appear peculiar often contain an 'elliptical' like morphology in WFC3. 
We show through several simulations that this larger fraction of elliptical-like galaxies 
is partially due to the courser PSF of WFC3, and that the 'elliptical' class very likely 
includes early- type disks. We also measure the merger history for our sample using 
CAS parameters, finding a redshift evolution increasing with redshift, and a peak 
merger fraction of ~ 30% at z 2 for the most massive galaxies with > 10^" , 
consistent with previous results from ACS and NICMOS. We also compare our results 
to semi-analytical model results and find a relatively good agreement between our mor- 
phological break-down and the predictions. We finally argue that galaxies classified 
visually as peculiar, elliptical and as peculiar ellipticals, all have similar structural and 
stellar population properties, suggesting that these galaxies are in a similar formation 
mode, likely driven by major mergers. 
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1 INTRODUCTION 

Explaining the origin of the Hubble sequence of galaxies has 
remained one of the most outstanding problems in extra- 
galactic astronomy. Massive galaxies in the nearby universe 
are largely in the form of ellipticals and spirals. These galax- 
ies tend to not be undergoing significant galaxy or star for- 
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mation processes, and are the end result of evolution and 
formation over the past ~ 13.7 Gyr. However, deep Hubble 
Space Telescope (HST) imaging using both optical and near- 
infrared (NIR) cameras, such as ACS and NICMOS, have 
shown convincingly that the galaxy population at higher 
redshifts, particularly at z > 2, is dominated by irregular 
and peculiar galaxies that often have no obvious similarity, 
in terms of structure, to lower redshift galaxies (e.g., Driver 
et al. 1995; Abraham et al. 1996; Conselice et al. 2005, 2008; 
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Lotz et al. 2004; Cassata et al. 2005, 2010; Overzier et al. 
2010; Cameron et al. 2010). In fact the structures and mor- 
phologies of galaxies, and how they evolve, are now recog- 
nised as one of the most important methods for understand- 
ing galaxies, as it reveals not only the empirical evolution of 
a new parameter, but furthermore correlates with physical 
processes, such as merging, occurring within galaxies (e.g., 
Conselice 2003; Conselice et al. 2003b; Cassata et al. 2005; 
Grogin et al. 2005; Trujillo et al. 2007). 

A major remaining issue in galaxy formation is solving 
the riddle of how the Hubble sequence came to be, and when 
this formation occurred. This problem of understanding the 
origin of the Hubble sequence can be defined most simply as 
uncovering when and how the disk and elliptical population 
in the nearby universe developed. This is a difficult problem 
since the population of disks and ellipticals is largely in place 
by 2 ~ 1 (Conselice et al. 2005) and to understand the 
formation of these systems one has to probe higher redshifts 
where the NIR is needed to resolve the location of most of 
the stellar mass in these galaxies. 

There are however some existing clues for how this for- 
mation occurs from deep Hubble Space Telescope imaging, 
particularly from observations with the NICMOS camera 
(e.g., Conselice et al. 2005; Conselice et al. 2011). These ob- 
servations show that, at rest-frame optical wavelengths, the 
types of galaxies seen changes gradually with time, such that 
disks and ellipticals are the dominate population at 2 < 1, 
while peculiars are the most common &t z > 2 (Conselice et 
al. 2005, 2008). The population of 'Hubble sequence' galax- 
ies roughly matches that of the peculiars sometime between 
2 = 1.5 — 2, which is also the peak of the star formation rate 
in the universe (e.g., Bouwens et al. 2009). 

There are at least two issues with these observations 
that still need to be addressed, and where HST Wide-Field 
Camera-3 (WFC3) observations will make a big impact in 
the next few years. One of these questions is understand- 
ing when the transformation between peculiar galaxies into 
disks and ellipticals occurs. The other major question is un- 
derstanding the physical mechanisms responsible for driving 
the formation of these Hubble sequence galaxies. For exam- 
ple, are ellipticals and spirals formed through the accretion 
of gas, such as in cold streams (e.g., Dekel et al. 2009), 
and/or through in-situ formation of stars within an exist- 
ing dark matter halo, and/or through mergers of existing 
galaxies (Hopkins et al. 2006)? 

Recent observations of the kinematic structures of high 
redshift galaxies show that there is a large diversity of types, 
but with many star forming galaxies having a large velocity 
dispersion compared to their rotation (Genzel et al. 2006). 
This is either an indication that these systems are large 
disks collapsing, or that they are undergoing some type of 
merging activity (Shapiro et al. 2008). Overall, kinematic 
studies of 2 > 2 galaxies find that up to 1/3 of these sys- 
tems have merger signatures in their kinematics, which is 
consistent with the fraction found from imaging (Conselice 
et al. 2003, 2008). However, these results are not entirely 
fool-proof, as other formation modes can mimic both the 
kinematic and morphological signatures. Furthermore the 
samples studied kinematically remain quite small, at a few 
10s of systems. The relative role of these various assembly 
processes at high-z remain an open and important question 




- / 

/ ... \ .... ^ ... . 
12 3 

Redshift (z) 



Figure 1. The rest- frame wavelength probed as a function of 
redshift for the two filters used in this study - Jno and Higo 
from WFC3 compared with the ACS filters of -B450, Vgoei *775 
and ^850- The vertical line at 2 = 1 denotes the lower Hmit to 
redshifts considered in this paper. The shaded region shows the 
region where the rest-frame optical is probed, showing that for all 
but near 2 ~ 1 one must use NIR imaging to study the rest-frame 
structures of these galaxies. 



regarding the mechanisms by which galaxy formation has 
occurred since z ~ 3. 

New deep NIR observations of the high-redshift universe 
can provide some answers to these questions by examining 
the first question of when the Hubble sequence was in place, 
and to begin addressing the formation mechanisms for galax- 
ies we see in today's universe. Observations with the WFC3 
of the Hubble Ultra Deep Field provides our first oppor- 
tunity to go beyond NICMOS observations that currently 
exist. We describe a preliminary investigation of this prob- 
lem using this early data and find that an apparently mostly 
smooth 'early-type' population appears to be the dominate 
one at high redshift, as compared with late- type or disk dom- 
inated galaxies. We further discuss the merger history and 
how it appears that high redshifts 1 < 2: < 3 were not a 
particularly hospitable time for the formation of prominent 
late-type disks. 

This paper is organised as follows. §2 gives an overview 
of the data and data sources we use, including previous work 
which we use to analyse our sample. §3 is a description of 
our structural measures, including CAS and Sersic profiles 
fits, as well as our measures of the spectral types of our 
galaxies. §4 discusses our results, including morphological 
fc-corrections, the types of galaxies found at high redshifts 
and a comparison to models, and how various classification 
methods compare with each other, while §5 is a summary 
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Figure 2. Stellar mass vs. redshift with the various galaxy types studied in this paper labelled. On the left we show the classifications for 
z > 1 galaxies carried out on ACS imaging of the UDF, and the right panel shows the classifications done using WFC3 at the rest-frame 
optical. The open blue circles are the galaxies classified visually as peculiars, the solid red boxes are ellipticals, SOs, and compacts, while 
the cyan triangles are ellipticals that appear to have a peculiarity. Disk galaxies are shown as green crosses for face on systems, with the 
edge-on systems displayed as a dot with a solid line. 



of the results of this paper. We use a standard cosmology of 
Ho — 70 km s^^ Mpc^^, and Qm = 1 — S7a ~ 0.3 throughout. 



2 DATA AND DATA SOURCES 

The primary data used in this paper is the WFC3 ob- 
servations of the Hubble Ultra Deep field located in the 
GOODS-South field (Giavalisco et al. 2004). We also use 
ACS and NICMOS imaging of the Hubble Ultra Deep Field 
(UDF) (Thompson et al. 2005; Beckwith et al. 2006) within 
this analysis, particularly when comparing with previous 
structural results (Conselice et al. 2008). The WFC3 field 
of view is 11 arcmin'^, with a drizzled pixel scale of 0.03 arc- 
sec per pixel, the same as for the ACS imaging. The filters 
we use in this study are the F160W (Hieo) and the FllOW 
(Hiio) with a depth similar to the optical down to roughly 
mAB ~ 29 using 16 orbits in the J-band and 28 orbits in H. 

The UDF ACS images are taken within the F435W 
(B435), F606W (Veoe), F775W (irrs), and F850L (zgis) 
bands. The central wavelengths of these filters, and their 
full- width at half-maximum are: F435W (4297, 1038 A), 
F606W (5907, 2342 A), F775W (7764, 1528 A), F850L 
(9445, 1229 A). In Figure 1 we show the rest-frame wave- 
length probed by the WFC3 filters we use, as well as those 
for the ACS, demonstrating the necessity of using the near- 
IR to probe the rest-frame optical at 2; > 1 (shown as the 
dashed area). 

The photometry and photometric redshifts we use are 
taken from our original study of the Hubble Ultra Deep 
Field (Conselice et al. 2008). The redshifts we use are pri- 



marily from Coe et al. (2006) who measure the photome- 
try of galaxies in the UDF within the BVizJH bands. The 
galaxies are detected with a modified version of SExtractor, 
called SExSeg. The photometry is PSF-corrected and aper- 
ture matched, removing the problem of matching magni- 
tudes at different wavelengths due to variations in the PSF. 
This high fidelity photometry is then used to derive photo- 
metric redshifts, and for the stellar masses we measure for 
our galaxy sample. 

The Coe et al. (2006) photometric redshifts are mea- 
sured using the photometric redshift techniques from Ben- 
itez (2000), based on the optical and NIR photometry from 
ACS and NICMOS. In addition to these photometric red- 
shifts we utilise 56 spectroscopic redshifts within the UDF 
field within our selection limits. All of the morphologies we 
determine in the near-infrared are examined with WFC3, 
and all the optical morphologies and structures with the 
ACS camera. 

As this is a morphological, and not a photometric study, 
we must have more stringent constraints on which galaxies 
we can examine. We limit our study to systems which have 
magnitudes zgso < 27, such that we are not biased by low 
signal to noise imaging. This final ACS based catalogue of 
2850 < 27 sources contains 1052 unique galaxies within the 
ACS UDF field. Of these 299 are found within the smaller 
field of view WFC3 imaging of the UDF within the redshift 
range 1 < z < 3. Note that although we use a zs5o limit 
to define our sample, we later limit our analysis in the var- 
ious analyses we carry out to a subsample of these. We do 
not include a discussion of lower or higher redshift galaxies 
in this paper. To run our morphological analysis requires 
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segmentation maps, and we ran SExtractor on the WFC3 
imaging to obtain these maps which we later use within the 
morphological analysis (see Conselice et al. 2008). 



3 STRUCTURAL MEASURES, STELLAR 
MASSES, AND SPECTRAL TYPES 

In this paper wc examine the morphological, structural and 
spectral properties of zsso < 27 galaxies in the UDF, as 
imaged in the WFC3 UDF pointing. We use an ACS «850 
selection as this was done in our previous analyses of the 
UDF (Conselice et al. 2008). This allows us to compare di- 
rectly how galaxy morphological and structural parameters 
vary with observed wavelength. Wc carry out several anal- 
yses of our WFC3 imaging using structural parameters, as 
well as stellar masses, and spectral types. We carry out these 
analysis in an independent way, thus lowering the effects of 
systematic errors on our overall results. 

We carry out our visual and CAS analysis in several 
steps to maximise the usefulness of the data, and to min- 
imise problems from contamination. We first create postage 
stamp images of each of our sample galaxies from our WFC3 
mosaic. These are created by cutting out a 10" x 10" box of 
the UDF surrounding each galexy, based on positions from 
the SExtractor catalog detections from Coe et al. (2006). 
Before this is done, the UDF WFC3 cutouts arc cleaned of 
projected nearby galaxies and stars through the use of the 
'segmentation map' produced by SExtractor (see Conselice 
et al. 2008). These segmentation maps are equivalent in size 
to the UDF image itself, with the difference being that it 
gives a numerical value for each pixel that reveals which 
galaxy it belongs to. These segmentation maps arc used for 
photometry, but they are also useful for removing nearby 
galaxies. The procedure we use is to replace pixels of galax- 
ies not being studied to the sky background with proper 
noise characteristics included. We then use these cleaned 
cutout images in our analysis. 

After examining our sample visually, we found that oc- 
casionally features remained near galaxies, and had to be 
manually removed by hand. There were also cases where 
large late-type galaxies with spiral arms brighter than their 
centres tended to be picked up by the program more than 
once, and these were manually noted when spotted. In the 
following sections wc describe our visual and quantitative 
analysis of these galaxy images within the UDF. 

Likewise, wc perform a parametric analysis on our data 
using the GALFIT package (Peng et al. 2002) to deter- 
mine parametric sizes and Sersic indices (§3.3). We use these 
Sorsic indices as another indicator of galaxy 'type' - i.e., 
what are likely progenitors of modern day disks and ellipti- 
cals. Furthermore, we also use the spectral energy distribu- 
tions of our galaxies to determine their spectral types, which 
we explain in §3.4. By using spectral types we have a third 
and largely independent way to determine galaxy types, or 
their progenitors, at higher redshifts. In the case, instead 
of structure, we are examining their stellar populations - 
namely, early- type, disk, and starburst/Im. 




WFC3 Class 

Figure 3. The relationship between the classification of the same 
galaxies within ACS and WFC3. The following key corresponds 
to the number schemes: = Elliptical, 0.5 = peculiar elliptical, 1 
= compact elliptical; 2 = early-disk, 3=late-disk , 6 = peculiar, 
8 = edge-on disk galaxy. Note that the points on this plot are 
randomly offset to show the relative number at each position. 



3.1 Visual Morphologies 

We use several different methods in this paper to exam- 
ine how the structures and morphologies of galaxies change 
with redshifts. First, one method involves a visual estimate 
of morphologies based on the appearance of our galaxies in 
the WFC3 imaging. The basic outline of how wc carry out 
our classification process is given in Conselice et al. (2005a), 
Conselice et al. (2007b) and ConseUce et al. (2008). The pro- 
cess is that we place each galaxy in the WFC3 imaging into 
one of nine categories: compact, elliptical, distorted ellipti- 
cal, lenticular (SO), early- type disk, late- type disk, edge-on 
disk, merger/peculiar, and unknown/too- faint. These clas- 
sifications are based only on how a galaxy looks, as seen in 
the observed Hieo-band within WFC3. These galaxies are all 
resolved within at least 3-4 times the effective radius. Wc do 
not use information such as colour, size, redshift, etc to de- 
termine our classification types. Note that with WFC3, the 
resolution is not generally very high for these distant galax- 
ies, and often galaxies look like 'blobs', and a classification 
type is sometimes more of an educated guess, and certainly 
prone to variation between classifiers. However, our classifi- 
cations were done consistently, as described below. First we 
give the main types in which our classifications were done. 

(i) Ellipticals : Ellipticals (Es) are centrally concentrated 
objects with no evidence for lower surface brightness, outer 
structures. Wc have 86 of these galaxies in our sample. In 
this paper, the meaning of 'elliptical' or 'early-type' hence- 
forth is likely not the same as a pure elliptical selection at 
low redshift. Simulations (§4.3) show that what we classify 
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Figure 4. The change in the CAS parameters as a function of redshift. We trace in this figure the change between the observed ~1.6 fim 
(Hieo) and ~0.9 fim (zgso) morphology. The symbol shape and colour denotes the galaxy type it denotes, as in Figure 2. The open blue 
circles are the galaxies classified visually as peculiars, the solid red squares are early-types, while the green crosses are disk galaxies. The 
horizontal line labels where there is no change between the measure structure in these two bands. Compared to morphological changes 
at 2 < 1 between rest-frame UV and optical, we find small changes at higher redshifts. 



as 'elliptical' is likely a mixture of true ellipticals and spiral 
galaxies dominated by a bulge and/or mergers in a more 
relaxed phase. 

(ii) Peculiar-Ellipticals : Peculiar ellipticals (Pec-Es) are 
galaxies that appear elliptical, but have some minor mor- 
phological peculiarity, such as offset isophotes, dual nuclei, 
or low-surface brightness asymmetries in their outer parts 
(58 systems). A full description of these galaxies is provided 
in Conselice et al. (2007). 

(iii) SOs: SOs are galaxies that appear to have a smooth 
disk with a bulge. These galaxies do not appear to have much 
star formation, and are selected in the same way nearby 
SOs are. Our sample has only a single system making its 
contribution very small. 

(iv) Compact - A galaxy is classified as compact if its 
structure is resolved, but still appears compact without any 
substructure. It is similar to the elliptical classification in 
that a system must appear very smooth and symmetric to be 
included. A compact galaxy differs from an elliptical in that 
it contains no features such as an extended light distribution 
or a hght envelope. The fact that these systems are common 
at 2 > 1 while ellipticals are more commonly seen at lower 
redshifts suggests that they are potentially drawn from the 
same population. (28 systems) 

(v) Early-type disks: If a galaxy contains a central con- 
centration with some evidence for lower surface brightness 
outer hght in the form of spiral arms or a disk, it is classified 
as an early- type disk. (14 systems) 

(vi) Late-type disks: Late-type disks are galaxies that ap- 
pear to have more outer low surface brightness disk light 
than inner concentrated light. (4 systems) 

(vii) Edge-on disks: disk systems seen edge-on, and whose 
face-on morphology cannot be determined, but is presum- 
ably an SO or spiral. (9 systems) 



(viii) Peculiar/irregular: Peculiars and irregulars are sys- 
tems that appear to be disturbed, or peculiar looking, in- 
cluding elongated/tailed sources. These galaxies are possi- 
bly in some phase of a merger (Conselice et al. 2003a), or 
are dominated by star formation (89 systems). 

(ix) Unknown/too- faint: If a galaxy is too faint for any 
reliable classification it was placed in this category. Often 
these galaxies appear as smudges without any structure. 
These could be disks or ellipticals, but their extreme faint- 
ness precludes a reliable classification. We find that there 
are no systems in this category using our selection, although 
unclassifiable systems are found in the UDF ACS imaging. 

Furthermore, we find that nine of our galaxies are not clas- 
sifiable as they are unresolved, appearing as point sources. 
This may be misidentified stars, QSOs, or a mixture of these. 
They are not included in further discussions in this paper. 
Furthermore, we performed these same classifications on the 
ACS Z850-band imaging of the UDF, in the same manner as 
for the Hi6o imaging, for comparison purposes. 



3.2 CAS Parameters 

In addition to visual morphologies, we utilise the CAS (con- 
centration, asymmetry, dumpiness) parameters to probe the 
structures of our galaxies quantitatively. The CAS parame- 
ters are a non-parametric method for measuring the forms 
of galaxies on resolved CCD images (e.g., Conselice et al. 
2000a; Bershady et al. 2000; Conselice et al. 2002; Conselice 
2003). The basic idea behind the CAS system (Conselice 
2003) is that galaxies have light distributions revealing their 
past and present formation modes. This system can also be 
used to find specific galaxy types, as defined in the nearby 
universe, such as ellipticals and mergers/peculiars. For ex- 
ample, the selection A > 0.35 locates systems which are 
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Figure 5. The relative distribution of galaxy types in the WFC3 Hubble Ultra Deep Field for systems selected with 2350 < 27 as 
measured within three different images. The left panel shows the distribution of our classifications using the ACS imaging of the UDF 
(Conselice et al. 2008). The middle panel shows our new WFC3 UDF classifications and the right panel shows these same classifications 
done on a WFC3 image convolved with the ACS point spread function. Labelled on each panel are disks, ellipticals/SOs, compact galaxies, 
and peculiars. Note that the compact galaxies become an important population at higher redshifts. 



nearly all major galaxy mergers in the nearby universe (e.g., 
Conselice et al. 2000b; Conselice 2003; Hernandez- Toledo et 
al. 2005; Conselice 2006b). 

The basic procedure for measuring the CAS parame- 
ters is essentially the same as that presented in Conselice et 
al. (2008). First the galaxy is cut out from the main UDF 
WFC3 image from which the radii and the CAS parame- 
ters are measured. We measure the Petrosian radii for each 
galaxy from our images. This is the radius defined as the lo- 
cation where the surface brightness at a given radius is 20% 
of the surface brightness within that radius (e.g., Bershady 
et al. 2000; Conselice 2003). The exact Petrosian radius we 
use to measure our parameters is 

iiPctr = 1.5 X r(?7 = 0.2), 

where r{r] = 0.2) is the radius where the surface brightness 
is 20% of the surface brightness within that radius. 

We use circular apertures to measure our Petrosian radii 
and quantitative parameter estimation. We begin our esti- 
mates of the galaxy centre for the radius measurement at 
the centroid of the galaxy's light distribution. Based on var- 
ious test and simulations, we have found that the resulting 
radii do not depend critically on the exact centre, although 
for other parameters this can be more important (Conselice 
et al. 2000; Lotz et al. 2004). For more detail see Bershady 
et al. (2000), Conselice et al. (2000), Conselice (2003) and 
Conselice et al. (2008). 

3.2.1 Asymmetry Parameter, A 

The basic method for measuring the asymmetry is by taking 
an original galaxy image and rotating it 180 degrees about 
its centre, and then subtracting these two images (Conselice 
1997). Important corrections are done for background hght, 
and the radius has to be measured carefully, as explained in 
detail in Conselice et al. (2000a, 2002). Details such as how 



centring is done is also described in Conselice et al. (2000a). 
The mathematical definition is 

. . /E|7o — /i8o|\ . / S|i3o — -B180I \ 

^ ^ \ nh\ ) ^ [ nio\ ) ' 

where Jo is the original image pixels, Jigo is the image after 
rotating by 180° . The background subtraction using light 
from a blank sky area, called Bo, are critical for this pro- 
cess, and must be minimised in the same way as the original 
galaxy itself. Lower values of A imply that a galaxy is largely 
symmetric, which tends to be found in early type galaxies. 
Higher values of A indicate an asymmetric light distribu- 
tion, which are usually found in spiral galaxies, or in the 
more extreme case, merger candidates. 



3.2.2 Light Concentration C 

The concentration index (C) is a measure of the intensity of 
light contained within a central region defined by a curve of 
growth radii, compared to a larger curve of growth radii. The 
exact definition is the log of the ratio of two circular radii 
which contain 20% and 80% (r2o, '"so) of the total galaxy 
flux, 

C = 5 X log . (2) 

This index is sometimes called C28. A higher value of C in- 
dicates that a larger amount of light in a galaxy is contained 
within a central region. Our measurement of the concentra- 
tion correlates well with the mass and halo properties of 
galaxies (e.g., Bershady et al. 2000; Conselice 2003). Nearby 
galaxy values for the concentration index are C = 2 — 3 for 
disks, C > 3.5 for massive ellipticals, while peculiars span 
the entire range (Conselice 2003). 
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3.2.3 dumpiness Parameter, S 

The dumpiness (5') is a parameter for describing the frac- 
tion of light in a galaxy contained in clumpy light concentra- 
tions. Galaxies that appear 'clumpy', such as star forming 
systems, have a relatively large amount of light at high spa- 
tial frequencies, and high S indices. Smoother galaxies, such 
as ellipticals contain smaller S values, dumpiness can be 
measured in a number of ways, the most common method 
used, as described in Conselice (2003) is. 



S = 10 X 



'S{Bx.y - B" 



(3) 



where, the original image Ix,y is blurred to produce a sec- 
ondary image, Ix,y This blurred image is then subtracted 
from the original image leaving a residual map, contain- 
ing only high frequency structures in the galaxy (Conselice 
2003). To quantify this, we normalise the summation of these 
residuals by the original galaxy's total light, and subtract 
from this the residual amount of sky after smoothing and 
subtracting it in the same way. The size of the smoothing 
kernel a is determined by the radius of the galaxy, and is 
a = 0.2-1.5 X r{ri = 0.2) (Conselice 2003). Note that the cen- 
tres of galaxies are removed when this procedure is carried 
out, typically the central l/8th of the Petrosian radius. 



3.3 Parametric Fitting — Sersic Indices 

While the non-parametric measures are one way to examine 
and quantify galaxy structure, utilising parametrised fits to 
light distributions of galaxies through an analytical form is 
another popular method for understanding the structural 
forms of galaxies. We fit to our WFC3 Higo-band imaging a 
single Sersic profile of the form: 



E(r) = Eeexp(— K[(r/re 



1]) 



(4) 



We use the GALFIT package to fit equation (4) above to de- 
termine the Sersic index (n) and the effective radius, (Rc). 
GALFIT was used to generate 2-D models of the light pro- 
files of our sample galaxies using a single Sersic function. 
The position, initial guess parameters for GALFIT (total 
mag. Re, ellipticity, position angle, etc.) were taken from the 
SExtractor catalog. The background value is kept fixed at 
the value determined by SExtractor in the area around each 
galaxy. The constraints used when performing the fitting in- 
cludes requiring that the centroid remains within 3 pixels of 
the input position, and the Sersic index, n, does not exceed 
a value of n = 10. The PSF for convolving the analytical 
model was generated from the stars in UDF field using the 
PSF tasks in IRAF/DAOPHOT package, and then normal- 
ising these. When there are neighboring objects surrounding 
a specific object of interest, we perform simultaneous fits to 
the neighbours, which has been found to work better than 
masking the objects based on some isophotal threshold. In 
this paper we mostly make use of the Sersic index, n, as 
a measure of the 'disk-like' or 'elliptical-like' nature of the 
galaxies within our sample. Typically galaxies with fitted 
n > 2 are elliptical like, and n < 2 are disk-like in a tra- 
ditional interpretation compared with nearby galaxies (e.g., 
Trujillo et al. 2007; Buitrago et al. 2008). With GALFIT 
we also measure the axis ratio of each system, giving an 
ellipticity of e = 1 — b/a. 
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Figure 6. The distribution of galaxy types within the ellipti- 
cal category in Figure 5 divided into normal ellipticals, compact 
ellipticals and peculiar ellipticals. 



While it remains to be determined whether the kine- 
matic interpretation of objects with rest-frame morphologies 
of a 'disk' and 'spheroid' nature match, the Sersic index is 
a good rough guide to the overall shape of the light distri- 
bution within galaxies at different redshifts. These fits were 
also done using WFC3 data in Cassata et al. (2010) to in- 
vestigate the size evolution of the most massive gala^xies at 
2: < 3, finding that they are indeed compact, as has been 
found in previous studies. 



3.4 Spectral Types 

Each galaxy in our sample has a morphological type which 
can be determined by the CAS parameters, by visual esti- 
mates of morphology, or by the measured Sersic index. We 
furthermore utilise in this paper another method of deter- 
mining galaxy types through the relative ages of their lu- 
minosity weighted stellar populations through fitting model 
spectra to observed spectral energy distributions (SEDs). 

These observed spectral energy distributions are taken 
from the wealth of ancillary photometric data in the UDF 
field. In summary we use TFITed magnitudes, including 
the ACS optical, Spitzer IRAC data, U-band ground based 
photometry and deep IR photometry from ground-based 
telescopes. In detail these spectral types are based on 
imaging data from the VLT/VIMOS (i7-band), HST/ACS 
imaging (F435W, F606W, F775W, and F850LP bands), 
VLT/ISAAC (J, H, and Ks bands), an Spitzer/IRAC bands 
at 3.6, 4.5, 5.8, and 8.0/im. The spectral types used in these 
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Figure 7. Example simulations of nearby galaxies which arc originally observed at z ~ imaged to how they would appear at z = 2.5 
within the Hubble Ultra Deep Field F160W band as imaged with WFC3. Tiiese nearby galaxies are classified as early-type spirals (Sa 
and Sb). The typical sizes of these galaxies are several kpc in effective radii, and are at a variety of distances (see Conselice et al. 2000a). 
As can be seen, many of these systems would be difficult to classify as early-type disks, or even as galaxies which contain a disk in some 
cases. 



fits include the SEDs based on tiie Coleman et al. (1980) 
types: E, Sbc, Scd, and Im based on spectra from nearby 
galaxies. In addition, starburst templates from Kinney et 
al. (1996) are included in these fits. We utilise the spectral 
types from these fits to compare with our morphological 
types. Further details of the fitting code, the templates used 
can be found in Dahlen et al. (2010). 



4 RESULTS 

In this section we analyse the morphological properties of 

our sample galaxies at z > 1. These galaxies constitute the 
brightest systems as seen in the Hubble Ultra Deep field and 
axe amongst the most massive galaxies at high rcdshifts. We 
examine our sample in several ways. One is to simply look at 
the distribution of galaxy types at « > 1 which we have clas- 
sified visually. The other is to examine how the morphologies 
of our systems differ between the WFC3 imaging and that 
seen in the ACS imaging of the same sample as examined in 
Conselice et al. (2008). We also examine how morphology, as 
judged visually, compares with spectral-type classifications 
from SEDs. Finally, we examine the structural parameters 
of our systems and compare what we calculate for the CAS 
values and Sersic fits for our galaxies and the spectral-types 
and visual morphological estimates for these systems. 



4.1 Morphological fc-corrections 

In Conselice et al. (2008) we investigate the morphologi- 
cal fc-corrcction for galaxies seen in the UDF ACS field, 
particularly those at z < 1. Wo carried this out to deter- 
mine how galaxies look different in a quantitative way at 
different wavelengths, and then used a quantitative version 
of this correction to put higher redshift galaxies at « > 1 
at a pscudo rest-frame optical wavelength from which their 
merger fraction was derived (Conselice et al. 2008). In this 
paper, we are now able to derive what is the actual morpho- 
logical fe-correction between the i?i6o-band and the zsso- 
band for « > 1 systems, or rather how rest-frame optical 
and rest-frame ultraviolet morphologies differ. 

First, to put this problem in context, the fundamental 
issue with using ACS data to examine the morphologies of 
z > 1 galaxies is shown in Figure 1 where the rest-frame 
wavelength probed by different ACS filters, and our WFC3 
filters, are shown. This demonstrates that to probe the rest- 
frame optical light of galaxies at 2 > 1.3 we must use the 
observed infrared, and that the F160W (Hieo) filter provides 
rest-frame optical coverage for galaxies between 1 < z < 
3. In the rest-frame optical, structural parameters such as 
the CAS parameters are largely similar (Taylor-Mager et 
al. 2007; Lanyon-Foster et al. 2011). Henceforth, later in 
this paper, we utilize mostly the Hieo imaging to probe the 
morphologies and structures of our galaxy sample. 

There are several ways to compare morphological k- 
corrections and how they evolve at 1 < a: < 3. One way 
to do this is to compare estimates of morphology in the ob- 
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Figure 8. Example simulations of nearby galaxies which are originally observed at 2; ~ imaged to how they would appear at z = 2.5 
within the Hubble Ultra Deep Field F160W band as imaged with WFC3. These nearby galaxies are elassified as late-type spirals (Sc and 
Sd). As can been seen, several systems are edge on in these simulations. The typical sizes of these galaxies are several kpc in effective 
radii, and are at a variety of distances (see Conselice et al. 2000a). 



served 2:350 ACS band with the WFC3 Hieo band. This com- 
parison is shown in Figure 3. In general, galaxies along the 
1:1 line have the same classification, and those which differ 
significantly, further away from this line, have increasingly 
different classifications in ACS and WFC3 

A few outstanding deviations should be noted on Fig- 
ure 3. The first is that for the peculiar class, as found in the 
ACS imaging, roughly half (78 out of 174) are also classified 
as peculiars in the WFC3 imaging. The remainder are disks 
(15 out of 174, or 9%), or early- types/compact galaxies (73 
out of 174, or 42%). This implies that either roughly half of 
the peculiars seen at 1 < a < 3 in the rest-frame UV are 
actually elliptical-like in the rest-frame optical, or that the 
larger PSF of the WFC3 imaging is creating galaxies that 
look smooth (see §4.2.1). We conclude that a significant frac- 
tion, but not all, of this change in type between wavelengths 
is due to the larger PSF of WFC3 compared with ACS §4.3. 

Otherwise, for those galaxies classified as early-type in 
the ACS image, we find that almost all, or 75 out of 83 sys- 
tems are also classified as an elliptical of some form in the 
WFC3 imaging. The remaining systems are mostly peculiars 
as seen by WFC3. The other slight difference in the classifi- 
cation between ACS and WFC3 is for systems classified as 
peculiar within the WFC3 imaging. For these systems, we 
find that 78 out of 89 systems (88%) are also classified as 
peculiar in the ACS imaging. The remainder are found to 
be early-types and two disks in the ACS imaging. 

We furthermore investigate the quantitative morpho- 
logical fc-corrections for 2; > 1 galaxies by comparing how 



the quantitative morphologies in the ACS zsso-band com- 
pares with the same morphologies in the ffioo-band. Figure 4 
shows the difference in the CAS parameters for our sample 
of galaxies from 2; ~ 1 — 3. What we find is generally similar 
to the trends seen for similar rest-frame wavelength mor- 
phological differences probed for z < 1 galaxies in Conselice 
et al. (2008) (see Tabic 1 and compare with Tabic 1 in Con- 
selice ct al. 2008). That is, the morphological fc-corrcction, 
particularly for galaxies near 2; ~ 1 in the rest-frame opti- 
cal and rest-frame UV are fairly similar to what is found at 
2 > 1. 



As expected galaxies are generally more asymmetric 
and clumpy when viewed at shorter wavelengths compared 
to the Hi6o-band, and surprisingly less concentrated for the 
ellipticals and peculiars. Table 1 lists the quantitative mor- 
phological fc-correction for our galaxies between 2 = 1 and 
2 = 3, using the same quantitative method as in Conselice 
et al. (2008). These values are the difference in the measured 
H160 morphology and the ACS 2350 structural features. 

We discuss below each of the CAS parameters and how 
their values differ between the Hieo-band and the 2850- 
band. First, what we find from our measurements is that 
the concentration values are most different between the 
Hiao and 2850 bands for the ellipticals, in the sense that 
the ellipticals are more concentrated in the z85o-band than 
in the Hieo-band. This is different from what is found for 
nearby galaxies (e.g., Taylor-Mager et al. 2007) where early 
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Figure 9. Plot showing tiie distribution of our WFC3 galaxy sample in the UDF as seen in the concentration-asymmetry plane of CAS 
space. The lines on each plot denote the region in which different galaxy types, as seen in the rest-frame optical are found in the nearby 
universe (e.g., Bershady et al. 2000; Conselice 2003). The redshift range for each panel is shown and goes from 1 < 2 < 3. We show both 
the visual and spectral type classifications within this figure. The main visual symbol types for the points are the same as in Figure 1, 
such that the open blue circles are the peculiars, the red boxes are the early-types, the cyan triangles are the peculiar ellipticals, the green 
lines are edge-on disks and the green crosses are disks/spirals. We also show the spectral type classifications for the same galaxies in the 
following way: if the galaxy has a red open box surrounding it, it is classified as a early- type spectrally, open blue boxes are starforming 
systems, Im/starbursts, and those that have a green centre are classified spectrally as disk like or Sbc/Scd. 



types become more concentrated at longer wavelengths. We 
also find this in the ACS UDF imaging at the lower red- 
shifts (2 = 0.25 — 0.75) where we find a positive value of 
AC /AX = 0.43 ± 1.5. However, we also find that early- 
type galaxies, as classified visually, are less concentrated in 
terms of their stellar mass distribution than in zgso-band 
light (Lanyon-Foster et al. 2011). This is likely due to cen- 
tral star formation dominating the morphologies of these 
systems, even if they are classified as early-type. 



We also find that the peculiar galaxies within 2 = 1 — 3 
have a similarly lower concentration as do ellipticals at 
longer wavelengths (Figure 2; Table 1). Spiral galaxies how- 
ever, have a concentration difference near zero, and on aver- 
age there is no morphological fc-correction for these systems, 
likely due to the fact that these spirals are dominated by star 
forming regions that are easily visible in both the rest-frame 
UV and rest-frame optical. 

In terms of the asymmetry parameter, we find a slight 
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Figure 10. The relation between the asymmetry (A) and dumpiness (S) for our sample of WFC3 UDF galaxies in the Hieo band. 
The solid line with the bend in it shows the relationship between these two parameters as found for nearby galaxies that are not 
involved in mergers (Conselice 2003). The straight line show the relationship between these two parameters found at similar rest- frame 
wavelengths by Lanyon-Foster et al. (2011). This relationship is such that non-merging galaxies with a higher dumpiness has a slightly 
higher asymmetry, both due to star formation. Merging galaxies, where the structure is distorted due to bulk asymmetries from a merger 
have a larger asymmetry for their dumpiness. The point symbols are the same as in Figure 1 and the redshift range for each individual 
panel is listed. Note that galaxies which are generally more clumpy than asymmetric are classified as disk/spirals visually. 



Table 1. The change in CAS parameters from the rest- frame op- 
tical to near-UV as function of redshift. These values are defined 
such that A = X{Hiqo) — -^(^gso) 
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-0.75±1.30 
0.08±0.54 
0.00±1.19 


Spirals 




z = 1 - 3 


|£(Mm- 
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II (Mm- 




0.01±1.20 
-0.34±0.48 
-0.09±0.23 



on average higher dumpiness values in the ^gso-band than 
in the iifieo-band. The ellipticals have on average the same 
dumpiness, which is often zero, within the iifieo-band and 
the zgso-band. 

Overall, what we find is that at z > 1 there is less of a 
morphological fc-correction than what is seen in the nearby 
universe, or at least at z < 1. This was originally noted in 
the lower resolution imaging of the Hubble Deep Field North 
when the NICMOS and WFC2 data were compared (e.g., 
Dickinson et al. 2000; Conselice et al. 2005). The meaning 
of this is fairly clear - the stellar populations that make up 
the hght originating in the rest-frame optical is dominated 
by younger stars, or at least stellar populations that contain 
enough young stars to still be visible in the ultraviolet. This 
implies that we are not likely to learn much beyond our cur- 
rent understanding of this aspect of galaxy morphology in 
larger surveys such as CANDELS - however, this does lead 
to the possibility of utilising the morphological fc-correction 
to find more evolved stellar population or ultra-dusty galax- 
ies at these redshifts. 



change with respect to wavelength, which is less pronounced 
than what was found in Conselice et al. (2008) for systems 
at z — 0.75 — 1.25, comparing rest-frame 0.3 fim to 0.5 /xm. 
As would be expected, the peculiar galaxies and the spiral 
galaxies are more asymmetric on average in the 285o-band 
than in the ffieo-band at 1 < z < 3 (Table 1). However, we 
find that the elliptical galaxies are slightly more asymmetric 
in the il/ieo-band on average, although the difi'erence is not 
large. 

The dumpiness index has a similar pattern as the asym- 
metry index in that the peculiar and spiral galaxies have 



4.2 Galaxy Structure at 2 > 1 

With WFC3 we are able to examine the rest-frame optical 
morphologies of gala^xies at 2 > 1 in detail using the WFC3 
Jiio and H160 imaging. We use multiple approaches for this, 
including visual estimates of morphological types, as well as 
examining the CAS parameters and Sersic fits to our galax- 
ies. 

The first observation we discuss is the distribution of 
galaxy visual structure with morphology and stellar mass, as 
shown in Figure 2. This Figure shows how our 2 > 1 sample 
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Figure 11. Plot showing the relationship between the concentration parameter C and the Sersic index, n for galaxies divided into 
redshift ranges 1 < 2 < 2 and 2 < z < 3. Show on this panel as well are the visual classification types, as in Figure 1, such that the open 
blue circles are peculiars, the red boxes are the early-types, the cyan triangles are the peculiar ellipticals, the green lines are edge-on 
disks and the green crosses are disks/spirals. 



is classified in both the ACS and WFC3 imaging of the UDF. 
The differences in these classifications can be visualised by 
comparing the two plots based on ACS and WFC3 imaging. 
A discussion of the visual apparent morphologies for these 
systems is included in Conselice et al. (2008). 

Figure 2 shows that there is a broad distribution of 
types over the entire stellar mass and redshift range of our 
sample of galaxies. The most obvious difference between the 
ACS and the WFC3 classifications is that there are fewer pe- 
culiar galaxies in the WFC3 identifications than are seen in 
the ACS images. The peculiar systems seen in the ACS imag- 
ing mostly become peculiar ellipticals and E/SO/Compact 
systems when viewed in the rest-frame optical with WFC3, 
as discussed in §4.1. 

We quantify this change by examining the distribution 
of morphological types as a function of redshift, as plotted 
in Figure 5. This shows the remarkable result that the pecu- 
liars which dominate the galaxy population within the ACS 
imaging, which probes the rest-frame UV ai z > 1, drop by 
a factor of 1.5-2 in relative distribution at z > 1 when im- 
aged within the WFC3 Hieo-band. A large fraction of this 
change in types goes into the elliptical and compact popu- 
lations. We do not see a large increase in the number of disk 
galaxies at high redshifts, however, we do see more of these 
systems in the WFC3 imaging (see §4.1). 

We further plot the breakdown within the elliptical clas- 
sification in Figure 6. Roughly half of the 'elliptical' popu- 
lation shown in Figure 5 is in the form of pure ellipticals, 
i.e., systems with no peculiarity. A further roughly 20% of 
the overall class of 'ellipticals' are compact and the remain- 
ing systems are peculiar ellipticals. There is no strong trend 
within the 'elliptical' class within these subtypes, although 



the peculiar ellipticals do appear to increase slightly at lower 
redshifts. 

If indeed the structures of the bulk of these galaxies are 
elliptical at high redshift this has important implications 
for galaxy formation studies. However, one major issue we 
have to deal with is the fact that the PSF of the WFC3 is 
larger than ACS and therefore it is important to quantify 
the effect of this larger WFC3 PSF on the classification of 
these galaxies. To address this we carry out several simula- 
tions to demonstrate that our results are robust. In §4.2.1 we 
convolve the PSF of WFC3 with the UDF ACS image and 
reclassifying the same galaxies. We also determine our clas- 
sification biases through simulating lower redshift galaxies 
to how they would appear at higher redshift in §4.3. 



4.3 Resolution and Depth Effect Biases 

4.3.1 Effects of the WFC3 PSF 

In all investigations of the structures and morphologies of 
distant galaxies a general problem is understanding how dis- 
tant galaxies appear, and how they would look if we could 
image them as we do nearby galaxies. We address this issue 
in some detail in §4.3.2. Before examining this problem in 
depth, we first ask the simple question - is what we see in 
terms of the difference in structure between our galaxy im- 
ages in the ACS and WFC3 due to the different PSFs, or a 
real effect? 

The difference between the ACS and WFC3 imaging 
morphologies seen in §4.2 is not simply just a matter of 
wavelength differences, but also PSF and level of depth. In 
principle, any of these three can create morphological dif- 
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Figure 12. Figure showing the distribution of GALFIT Sersic indices for our visually classified galaxies within the WFC3 UDF. Show 
are two panels divided by redshifts, with the left hand panel galaxies in our sample found at 1 < z < 2, while the right panel shows 
galaxies at 2 < 2 < 3. For each visually determined galaxy type: elliptical (E), peculiar elliptical (pE), spiral/disk (Sp), cdgc-on disk 
(EO) and peculiar (pec), we show the distribution of Sersic indices. The vertical lines at n = 2 demonstrates the separation between 
galaxies that are 'disk-like' with n < 2 and those that are 'elliptical-like' with n > 2 (see Ravindranath et al. 2006; Buitrago et al. 2008). 



ferenccs (Consclicc ct al. 2008), and while the first one due 
to wavelength is a 'real' (i.e., not instrument driven) effect, 
if there are differences due to PSF or depth, tlien this is 
simply due to the way the data was taken. To address this 
issue, we convolve the ACS image of the UDF in the zgso- 
band where our original work was carried out (Consclicc et 
al. 2008) with the WFC3 PSF to determine how the mor- 
phologies of galaxies changes solely due to the slighter worse 
resolution of WFC3. 

To carry this out we reclassify our galaxies seen in the 
new convolved image in the same way we did on the orig- 
inal ACS and WFC3 images (Figure 5). We carry out this 
classification in the same way we did our original classifica- 
tions. This allows us to (ictcrininc, and separate, the contri- 
bution to the change in morphology as resulting from differ- 
ent PSFs between the ACS and WFC3 camera as opposed to 
differences produced from probing different rest-frame wave- 
lengths. 

What we find is that the disk galaxy fraction is similar 

between the ACS, WFC3, and the ACS convolved with the 
WFC3 PSF image, classifications. However, there is clearly a 
difference between the fractions of peculiars and early-types 
when comparing the ACS WFC3 PSF convolved image with 
the WFC3 one. While we find a large fraction of E/Compact 
systems within the WFC3 images at all redshifts, we do not 
see the same results in the ACS WFC3 convolved image. 
In fact, while the number of early- type systems in the con- 
volved ACS imaging increases from roughly a 30% fraction 
to about 45% after convolving with the WFC3 PSF. The 
peculiar fraction also drops from roughly 60% to 50%. This 
is significantly different from the WFC3 imaging classifica- 
tions where ellipticals make up 60% of the population and 
peculiars are at ~ 30%. This demonstration shows that over 



half of our classifications in the WFC3 imaging are either 
due to identifying the correct galaxy type, or are due to 
morphological fc-corrections. 

4-3.2 Galaxy Structure Simulations 

In the previous section wc examine how galaxies in the UDF 
ACS image would appear after convolved with the WFC3 
PSF. This approach is useful when comparing the WFC3 
and ACS classifications, but they do not account entirely 
for the effects of redshift on the appearance of galaxies as 
they change through redshift. To account for this, wc carry 
out a series of full simulations of taking galaxies at « ^ 
to how they would appear at higher redshifts. Within this 
discussion, we use simulations of nearby galaxies placed at 
a = 2.5 as a fiducial amount of change for nearby galaxies 
to within the redshifts of this study (1 < 2 < 3). 

These simulations axe based on the method of simulat- 
ing nearby galaxies at higher redshift as outlined in Con- 
selice et al. (2003). Similar simulations have been carried 
out when examining morphologies in ACS imaging in Con- 
sehce (2003) and Conselice et al. (2008), and for NICMOS 
and WFC2 imaging in Conselice et al. (2003) . These simula- 
tions are also used to determine how quantitative structural 
parameters change with redshift as outlined in Conselice et 
al. (2003, 2008). 

These simulations are based on real images, taken in op- 
tical B-band, for nearby galaxies, using the sample from Frei 
et al. (1996). This sample is discussed in terms of morphol- 
ogy in Consclicc et al. (2000a) and Bcrshady et al. (2000). 
Suffice it to say that this sample consists of bright nearby 
galaxies, mostly ellipticals and spirals. We do not consider 
merger simulations in this paper, although they have been 
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Figure 13. Figure showing the distribution of spectral types for our visually classified galaxies within the UDF. Show are two panels 
divided by redshifts, with the left hand panel galaxies in our sample found at 1 < 2 < 2, while the right panel shows galaxies at 2 < z < 3. 
For each visually determined galaxy type: elliptical (E), peculiar elliptical (pE), spiral/disk (Sp), edge-on disk (EO) and peculiar (pec), 
we show the distribution of spectral types for each. These types are such that: 1 = E, 2=Sbc, 3=Scd, 4=Im and 6=Starburst. The 
vertical lines show the division between early types at < 2, mid-types or disk at 2 — 5, and star forming galaxies, such as irregulars at 
> 5. 



discussed in detail in Conselice (2003), and will be the focus 
of future papers on this topic. 

The simulations arc done essentially by taking these 
nearby galaxies, £is imaged with the Lowell 1.5 m telescope 
and placing them as they would appear in our WFC3 imag- 
ing of the UDF. All aspects of the Hubble Space Telescope, 
and the WFC3 imager, are considered in this simulation, in- 
cluding its: aperture, baud width of different filters, through- 
put, pixel scale, read-noise, dark current, as well as using the 
empirically derived sky background (e.g., Giavalisco et al. 
1996; Conselice 2003). Furthermore, we also experimented 
with two PSFs - one calculated from Tiny Tim and another 
based on stacking real stars from WFC3 imaging in other 
fields, such as the GOODS Early Release Observations. 

The process in carrying out these simulations is outlined 
in detail in Conselice (2003) and we only give a summary 
of this here. The nearby galaxy image is first reduced in 
angular size to how it would appear at the simulated red- 
shift. Next the flux is reduced by (1 -I- z)* to account for 
surface brightness dimming. A background is then created 
which include all noise contributions, the major one being 
from the WFC3 background, which is at a level of ~ 0.5 e~ 
s~^ pix~^, which is around a factor of ten higher than for 
ACS. The galaxy is then placed into this simulated back- 
ground, and the WFC3 PSF is applied to the image. We 
increase the surface brightness of each image by 1 magni- 
tude, as it is known that distant galaxies are brighter than 
nearby galaxies by roughly this amount, with passive evo- 
lution with an exponentially declining star formation rate 
producing at least this much evolution. 

Figure 7-8 shows examples of these simulations, for 
galaxies both before and after the simulation. We only show 



in this paper the disk galaxies, as mergers and ellipticals 
nearly always are identifiable as such when simulated at 
higher redshift. We divide these disks in early-types (Sa 
through Sb; Figure 7) and late-types (Sc and Sd; Figure 8). 
The results are at times worrying for the simple reason that 
many of these galaxies would appear to be classified visually 
as ellipticals. In fact, the problem is more acute for the early 
type disks where the bulge component tends to dominate the 
appearance of the galaxy after it has been simulated. View- 
ing these galaxies as we did with ds9 it would be difiicult to 
always identify these as disks, especially as morphological 
classification is a subjective process, and different classifiers 
will classify the same galaxy in a different way. 

In general for the late-type (Figure 8) and to a lesser 
extent the early-type disks (Figure 7) the simulated galaxies 
appear to have a morphology that is mostly bulge. However, 
these systems often have fainter outer 'halo'-like material 
which can be interpreted, if identified correctly, as a disk 
system. However, it is possible that these systems could be 
identified as ellipticals with an outer envelope. In this paper, 
we were therefore careful to carry out our classifications of 
disks based on whether or not a roughly symmetrical outer 
portion of light was detected in an otherwise bulge system. 
Regardless, it is obvious that even with deep UDF WFC3 
imaging that much of the structure in spiral galaxies is lost. 
Clearly spiral arms can no longer be identified, and often the 
disk itself is much less prominent than the bulge, undoubt- 
edly this is due to the lower surface brightness of disks. This 
will be a major issue when interpreting shallower CANDELS 
data. 

We classify our simulated sample as we did our orig- 
inal, and find that within these conditions we would have 
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misclassified l/4th of the disk galaxies in the Frei samples 
as 'early- types'. This is similar to the fraction of galaxies we 
misclassify in the ACS imaging due to the larger PSF. The 
reasons for these is likely similar, and due to the larger PSF 
of WFC3. Wc account for this in our analysis of the overall 
evolution of types in §4.4. 

A major caveat to these simulations is that it is likely 
that disk galaxies in the early universe are potentially differ- 
ent from modern disks. We know that galaxies, even those 
with a disk-like morphologies and Sorsic indices indicating 
disks, are smaller by a factor of a few at 2; ^ 2.5 for at least 
the most massive systems (Buitrago et al. 2008). This will 
make it harder to identify disks if they arc indeed smaller 
and therefore more difficult to resolve. On the other hand, 
we know that the quantitative morphological distribution 
for massive galaxies is quite different at these redshifts than 
today, and that any disks at these redshifts will likely not 
have as prominent bulges as contemporary disks. Further- 
more, disks at high redshifts are potentially undergoing sig- 
nificant amounts of star formation, making their disks easier 
to identify and study. All of these features will have to be 
examined in more detail in future papers. 

4.4 Overall Results of Apparent Morphological 
Evolution 

In summary, the correction for PSF effects within WFC3 is 
at a level of 30-40% at most, as described in §4.3.1. This 
is verified by our analysis of simulations of nearby galaxies 
placed at high redshift (§4.3.2). In these simulations, we still 
find that 'carly-typc' morphologies are the dominate one at 
1 < z < 3. At most a correction would bring down the 
E/compact line in Figure 5 to roughly 40-50% - which is still 
much higher than what is implied with the ACS imaging, as 
well as different from previous results (Conselice et al. 2008). 
Therefore one major conclusion fi'oiii this work is that high 
redshift galaxies appear to be more early-type than disk- 
like. However, this is only for our mass range, which is from 
roughly log M, = 9 to just above log M, — 10. This trend 
is potentially different for galaxies with log M* > 11 as 
seen in the GOODS NICMOS Survey (Buitrago, Trujillo & 
Conselice 2011). 

However, as we show in the remainder of this paper, 
these galaxies, while appearing to be early-types are any- 
thing but analogs of early-type, in terms of their stellar 
populations, and they are likely undergoing significant evo- 
lution. Already, this can be seen by the fact that there is a 
morphological fc-correction for WFC3 identified early- types, 
which appear as distorted systems in the rcst-framc UV as 
probed by the ACS camera. We discuss the structures of 
these early-types as well as their spectral types, demon- 
strating the ongoing evolution of the Hubble sequence at 
l<z<3. 

4.5 Quantitative Structure and Spectral Types 

4-5.1 CAS and Sersic Classifications 

In Figure 9 wc show the locations of our sample in 
concentration-asymmetry space and label our galaxies ac- 
cording to both their spectral and visual classification type. 
This plot is best understood by looking for boxes that agree. 



or otherwise, in terms of the two colours plotted at each 
point. A blue box is for a starburst/Im spectral type, and 
a red box is for an E spectral type, while a green central 
dot is for Sbc/Scd spectral types. If the two colours are the 
same then the two ways of classifying that galaxy through 
visual and SED spectral-types are the same. If they differ 
then they do not agree. 

The most remarkable feature of this diagram is that 
there is a significant amount of systems which are classi- 
fied as having different spectral and morphological types 
throughout. What can be seen however, is that there is some 
relation between spectral galaxy type and position on the 
concentration-asymmetry plane. The starburst/Im classifi- 
cations are nearly all in the higher asymmetry part of the 
distribution at z < 2. This however, changes at 2; > 2 where 
the starburst/Im types arc found throughout the concentra- 
tion/asymmetry plane. This is another indication that the 
quantitative measurements of structure at 2: > 2, particu- 
larly for lower asymmetry systems, cannot be interpreted in 
exactly the same way as a lower redshifts (e.g., Conselice et 
al. 2008). 

Furthermore, as Figure 9 shows, there is not always a 
good agreement between visual estimates of morphology and 
their positions in the CA space. The disk galaxies, both 
those classified as such visually and through spectral en- 
ergy fits to SEDs are often different objects. Furthermore 
the peculiar galaxies have a high asymmetry at each red- 
shift, and often the peculiars are the majority of the sys- 
tems in the merger region of CAS space, with A > 0.35. 
However, those galaxies classified visually as early-types arc 
often found in the region of CA space occupied by the lo- 
cation of late-type galaxies in the nearby universe. This 
does not imply that these systems are late- types, but that 
they have asymmetry values that are larger, with around 
< A >~ 0.2. This is significantly higher than the average 
asymmetry for early types in the nearby universe which have 
values < A >= 0.02 ± 0.02, which is more than 10 a differ- 
ent from higher redshift galaxies. This clearly shows that 
the visually identified ellipticals are in a formation mode. 

We furthermore show on Figure 10, the asymmetry- 
dumpiness plane of the CAS space. Many objects are too 
smoothed by the WFC3 PSF to measure this parameter (sec 
Conselice 2003). However, those systems which are the most 
clumpy, as compared to their asymmetry, are the galaxies 
identified visually as disk systems. This is due to the fact 
that these visually identified disks are undergoing intense 
star formation in clumps, and potentially, using a compari- 
son of asymmetry vs. dumpiness is a powerful way to iden- 
tify these systems in WFC3 data. 

We can get a further idea regarding the distribution 
of the properties of various high redshift Hubble types by 
examining the Sersic fitted n index vs. the concentration 
parameter for our systems as labeled by spectral type. Fig- 
ure 11 shows this, where we plot the spectral type as different 
symbols. In general, we find that galaxies with higher Sersic 
indices have higher concentration values, with most systems 
at n > 2 having a concentration value C > 2.7, typically 
the value for early-types in local universe (Conselice 2003). 
However, clearly there is significant scatter in these values, 
and the galaxies identified visually as early-types (red boxes) 
are in every location of this diagram. The disk galaxy sys- 
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tems (the green crosses) however show that at least the disk 
galaxy systems are found at lower Scrsic indices of n < 2. 

Figure 12 shows the distribution of the Sersic indices 
for our sample, showing the range of values for each visually 
classified type, and that most visually classified ellipticals 
have lower Sersic indices. This demonstrates that at least 
a fraction of the early-types, or the progenitors of early- 
types, have the expected Sersic morphology, but still have 
enough residual star formation to have spectral types such 
as Scd or starburst (Figure 13; §4.5.2). Figure 14 on the 
right hand side, furthermore shows that much of the n > 2 
systems have a starburst spectral-type, and thus we may be 
seeing the epoch of the formation of the bulk of early types 
at z > 2 when finding bona-fidc early types as a non-star 
forming, morphologically identifiable early-type, will be rare 
at these times (Bauer et al. 2011). 

4-5.2 Spectral Type Classifications 

When examining galaxies and classifying them into different 
types, there are several ways to proceed. The traditional way 
is to examine systems in optical light and to classify the 
galax;y onto the Hubble sequence or a variation of the Hubble 
sequence. For nearby galaxies there is a good correlation 
between the Hubble type and physical parameters of the 
galax;y, such as the colour, star formation history, size, etc. 
(e.g., Roberts & Haynes 1994; Conselice 2006a). However, 
this relationship may not hold at higher redshifts. Therefore, 
it is desirable to classify gala:xies into spectral types, as well 
as visual estimates of morphology. 

Spectral types can be obtained through high S /N spec- 
tra, or through fitting of the spectral energy distributions 
of the galaxies. The approach that we use in this paper is 
to fit spectral energy distributions to the broad-band pho- 
tometry of our sources and to derive which is the best fitting 
spectral type - elliptical, early-type spirals, late-type spirals, 
irregulars and starbursts (§3.4) . This is similar to our mor- 
phological scheme in principle, although there is no reason 
why these spectral-types should match morphological ones 
at 1 < 2 < 3 as well as they do at lower redshifts. 

First, we give some rmmbers on the overlap of these 
types. At 1 < 2 < 2, we find that for systems identified as 
disk galaxies and edge-on disk galaxies visually, 86% (83%) 
of systems have a spectral- type consistent with a disk, i.e., 
Sbc/Scd types (where the second number is parenthesis is 
the edge-on disk fraction). For the remainder, we find that 
7% (17%) of disks have an early-type spectral-type, and 7% 
(0%) have spectral types that are star forming/starbursts. 
What this shows is that our visual estimates are roughly 
correct in identifying which systems have a mix of stellar 
populations - i.e., a mixture of an old and a younger compo- 
nent. At higher redshifts (2 < « < 3) we find for both types 
of disks that 100% are identified as Sbc/Scd spectral types. 

For systems identified as visual early-types in WFC3 
imaging, we find a more mixed outcome. At 1 < z < 2 wc 
find that only 10% of these systems would be identified as 
E-spectral types. That is, only a very small fraction of the 
WFC3 imaged UDF galaxies we find visually as early-types 
have a spectral shape suggesting that they contain evolved 
stellar populations. The bulk of the visually classified early- 
types are found to be consistent with a Sbc/Scd spectra type 
(56%) at 1 < z < 2, suggesting a mixture of old-|-young stel- 



lar populations, or the fading of an older stellar population 
that did not form in one complete burst. Interestingly, 34% 
of the visually classified early-types at 1 < « < 2 are consis- 
tent with being starbursts or Im systems - that is galaxies 
which are dominated by star formation. This rises to 62% 
starburst spectral-types for visually classified early-types at 
2 < z < 3. We will return to a discussion of this issue in §5 
and what it implies for early-type formation. 

The two final visual categories we discuss are the pe- 
culiars and the peculiar early-types, which are systems that 
resemble early-types in some formation state and display a 
visual peculiarity. What we find is that at the higher red- 
shift range, 2 < z < 3, these two types have similar spec- 
tral energy distributions - the visually classified peculiars 
are 62% starburst/Im spectral type, and 38% Sbc/Scd type. 
Very interestingly, this is just slightly less than the fraction 
of visually classified early-types which are star forming. As 
for the peculiar early-types, we find a very similar fraction 
of 60% Im/starburst spectral- type and 40% Sbc/Scd. Based 
on this, it appears that the peculiar/E/pE types may have 
a similar origin, and this is reflected in the similar spectral 
types of all these types. 

At lower redshifts, 1 < z < 2, we find that the fraction 
of peculiars which remain starburst/Im systems drops to 
35%, with 10% having an E spectral type, and 55% having 
a Sbc/Scd spectral type. This is very similar to the break- 
down for the ellipticals. This shows that peculiar galaxies 
arc not all starbursting systems, but likely are distorted due 
to perhaps the merging of existing galaxies which contain a 
mixture of old and young stellar populations. The peculiar 
ellipticals at 1 < « < 2 are 53% starburst/Im spectral types 
and 47% Sbc/Scd spectral types, showing a similar pattern 
to the peculiars, but with even more systems having a star- 
burst origin. 



4.6 Tracing Galaxy Formation and Actual Types 

There are two goals to studying galaxies in terms of their 
structure in this paper. These are understanding when and 
how the Hubble sequence was put into place. Clearly, as we 
have just shown, the Hubble sequence, as we know it today, 
does not exist at higher redshifts, particularly at z > 2. The 
question becomes - how can we identify galaxies at higher 
redshifts that match the properties of nearby Hubble types? 
This goes beyond morphology and asks the question - when 
do ellipticals, in the sense of smooth, nearly structureless 
galaxies dominated by old stellar populations, first form? 
Related to this is the formation of galaxies with disks. The 
former problem of the ellipticals is easier to probe, in princi- 
ple, as finding a galaxy with a smooth morphology and older 
stellar populations is more straightforward than determin- 
ing if a galaxy is a rotating disk. Wc also investigate whether 
we can use structure to determine the formation modes of 
galaxies, where we consider the merger history in §4.6.3. 

4.6.1 A 'Pure' Elliptical Evolution 

First, we discuss the fraction of pure ellipticals within the 
UDF. We define pure ellipticals as galaxies which have a 
spectral-type consistent with an early-type, and a morphol- 
ogy consistent with being a modern elliptical, which we se- 
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Figure 14. Plot showing the relationship between the concentration parameter C and the Sersic index, n for galaxies divided into 
redshift ranges I < z < 2 and 2 < 2 < 3. Show on this panel as well are the spectral classification types, with the following labelling 
scheme: the solid red open boxes are classified as early-type spectrally, open blue circles are starforming systems, the Im/starbursts, and 
the green crosses are those classified spectrally as disk like, or Sbc/Scd. 
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Figure 15. The merger fraction as a function of redshift and stellar mass for our WFC3 sample (blue triangles) compared with data 
from the Hubble-Deep Field North NICMOS (red squares) and the ACS Ultra-Deep Field (black circles). Shown in this figure are the 
merger fractions at mass limits of 8 < log M« < 9, 9 < log M« < 10, and log M» > 10. The green crosses show the evolution of the 
merger fraction for a combined UDF and HDF-N sample. 



lect asA<0.2. At2 = the elliptical fraction for galax- 
ies with stellar masses M, > 10^° Mq is > 20% (Conselice 
2006a). With our selection, we find a much smaller fraction 
at higher redshifts. 

When we examine the evolution of this pure elliptical 
fraction with redshift we find that there are no pure early- 
types ai z > 2. We find a total of six systems at 1 < 2; < 2 



giving a fraction of 3.2±2.3% at 1.5 < z < 2.0 and 4±2% at 
1 < z < 1.5. We find a higher fraction if we limit our selec- 
tion to M, > 10^° Mq systems = 14±14% at 1.5 < z < 2 
and 23±13% at 1 < z < 1.5. This demonstrates that mas- 
sive galaxies are more likely to be elliptical than lower mass 
ones, which is also seen in the nearby universe (e.g., Con- 
selice 2006a). This also shows that at least within the very 
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small field of view of the UDF, there are no pure ellipti- 
cal galaxies at z > 2. However, most certainly progenitors 
of pure ellipticals exist at these redshifts, and the fact that 
there are massive galaxies, even within our own sample, at 
z > 2 shows that these galaxies arc in an active assembly 
mode, and can be identified and studied. Future studies of 
larger fields, especially within the CANDELS programme, 
will address these issues, and put better constraints on the 
number of passive ellipticals at 2 > 2 

4-6.2 The Disk Galaxies 

Tracing disk galaxies in a similar way as the ellipticals is not 
as straightforward. However, the disk galaxies at the magni- 
tude and stellar mass limits we use are extremely common 

in the nearby universe, with a fraction of 60-75%. Tims, just 
from this, a significant number of the distant galaxies we see 
must be progenitors, in some form, of disk galaxies in today's 
universe. The issue with disks is that the spectral-type is a 
mixed population of older and younger stars, which is not 
an unusual situation at high redshift, and does not necessar- 
ily imply that these systems are always disks - consisting of 
an older bulge with a younger disk. Structurally disks can 
be identified by having a modest asymmetry and concen- 
tration. When we examine systems within our sample with 
C < 3 and 0.2 < ^4 < 0.35, and spectral types consistent 
with nearby disks, we find a roughly even fraction of disks 
of 7-15% between redshifts of 1 < a < 3, which is similar to 
the visual fraction (Figure 5). 

However, when we use these limits we find that the 
galaxies we select are not always identifiable as disks vi- 
sually, in fact they often are not. We in fact classify most 
of these systems as peculiars or compact galaxies, showing 
that the apparent morphology of galaxies can be quite dif- 
ferent from what their structure and stellar populations are 
revealing. 

An alternative method for understanding the limits of 
the number of disks we have is by examining the fraction of 
systems with large inclinations, which arc likely disks seen 
edge-on. Through the use of our measured ellipticities with 
GALFIT, we find that roughly 30% of our selected systems 
have an axis ratio b/a < 0.3. If all of these are edge on sys- 
tems, then every galaxy at these redshifts from z — 1 — 3 
would have to be a 'disk' galaxies. While this is obviously 
unlikely the case, it shows that we are likely misidentifying 
some disks, although many of these disks are likely in some 
formation mode and cannot be identified as such. This is- 
sue was furthermore examined recently for massive compact 
quiescent galaxies at a ~ 2 by van der Wei et al. (2011) who 
found that 65±15% of massive compact systems are likely 
disks in some form. While it is unlikely that every galax- 
ies in our sample is a disk in some form, it docs show that 
we are very likely under representing the number of disk 
galaxies in our sample by using either visual morphologies 
or parametric parameters (see also Buitrago et al. 2011, in 
prep). 

4.6.3 The Merger History 

In this section we revisit the measurement of the merger 
history which was used in Conselice et al. (2008) with ACS 



UDF data, and Conselice et al. (2003) using NICMOS imag- 
ing of the HDF-N to trace the merger history for galaxies 
as a function of stellar mass at ^; > 1. Both of these studies 
found similar results, that the merger history is highest for 
the most massive galaxies at z ~ 2.5, and declines steeply at 
lower redshifts, while lower mass systems have a more grad- 
ual drop off in their merger fraction. Conselice et al. (2008) 
and Conselice (2006b) investigate, using these merger frac- 
tions, the likely merger rate, and the contribution of merging 
to building up the stellar masses of galaxies over time. 

We use WFC3 data from this paper to make a new 
measurement of the merger history with WFC3 data at ^; > 
1, and compare this with previous work. The errors on these 
measurements are large simply because of the small number 
of galaxies used in this study. Future work using data such 
as CANDELS will address the merger history issue in much 
more detail. For now however, we simply describe the merger 
history, and how it agrees with previous results. 

First we describe our measured merger fractions based 
on the CAS parameters. This is a topic with a large history 
and background, but the most relevant papers for measuring 
the merger history with CAS are Conselice (2003), Conselice 
ct al. (2003), Conselice et al. (2008) and Conschcc, Yang & 
Bluck (2009). As described in Conselice et al. (2000a,b) and 
Conselice (2003) one method for finding mergers is to use 
the conditions in rest-frame optical light (cf. Taylor et al. 
2007; Lanyon-Foster et al. 2011 for UV), 

A>0.35&^>S', (5) 

such that the asymmetry is higher than A = 0.35, and the 
asymmetry value is higher than the dumpiness value, A > S. 
This ensures not only that a galaxy is roughly 1/3 asymmet- 
ric in terms of its light distribution, but also that this light 
does not come from clumpy light, such as star forming re- 
gions. This method works well at lower redshifts (Conselice 
2003), yet with the larger star forming clumps seen within 
high-z galaxies (e.g., Elmegreen et al. 2005) it remains pos- 
sible that some galaxies with large star forming complexes 
are being misidentified as a merger. 

We show the merger fraction using our new WFC3 UDF 
imaging in Figure 15, finding a very similar merger fraction 
as before using ACS. One way to quantify this merger his- 
tory is through various fits of the merger fraction history. 
The traditional power-law format (Conselice et al. 2003a,b; 
Bridge et al. 2007; Conselice et al. 2009; Bluck et al. 2009) 
is given by: 

/m(;?) = /o X (1 + zr (6) 

where fm{z) is the merger fraction at a given redshift, /o is 
the merger fraction at 2; = 0, and m is the power-law index 
for characterising the merger fraction evolution. An alterna- 
tive way to characterise the merger fraction evolution is the 
Press-Schechter formalism (Conselice et al. 2008) - which is 
essentially a combined power-law-|-exponcntial form. How- 
ever, we do not consider this form in this paper, but will in 
future investigations of the merger history. 

By fitting equation 6, we find a similar parameterisa- 
tion of the merger history as we found for galaxies with an 
m exponent of m ~ 5.5, for massive galaxies and rn ~ 2.2 
for lower mass systems. We also show on Figure 15 the best 
fitting exponential/power-law form for the merger history. 
Integrating the merger rate, using a time-scale of 0.4 Gyr 
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(Conselice et al. 2008) we find that on average there are 
4.3I0 g major mergers at 2: < 3 for the most massive galax- 
ies with M, > 10^^ Mq . This is enough merging to roughly 
double the mass of these galaxies (Conselice 2006b). In the 
future, we will utilising more wide-field imaging from CAN- 
DELS and other WFC3 surveys to obtain a more accurate 
and definite merger history for these galaxies. 

4.7 Comparison to Models 

In this section we carry out a comparison between our re- 
sults, and models of galaxy formation which can predict the 
evolution of galaxy morphology as a function of redshift. The 
models we compare with are the Galacticus semi-analytical 
models from Benson & Devereux (2010) and Benson (2010) 
which do not resolve the morphologies of galaxies as some 
numerical models do. However, these models follow the for- 
mation of galaxies in a cosmological context, and can be 
used to predict how spirals, ellipticals and galaxy mergers 
evolve through time based on the physics which is thought 
to form these systems^ 

Currently, the idea is that galaxies are formed from 
mergers and accretion. Bulges and ellipticals are formed in 
merger events that have dynamically cooled and relaxed, 
while disks and spirals are formed by the smooth accretion 
of intergalactic gas which is then converted into stars over 
some time period. These models predict these quantities, as 
well as the time since the last major merger for each galaxy 
at a given redshift. 

We use these models to predict the morphological evo- 
lution of galaxies in a cosmological context, and therefore 
whose formation is influenced by the relative distribution of 
dark matter, dark energy and the temperature of the dark 
matter. These models are such that a measure of the bulge 
to total ratio can be derived based on the history of spheroid 
and disk assembly. Thus, while these models do not resolve 
galaxy morphology, they can predict what the likely mor- 
phological distribution is based on the physics of the galaxy 
formation. 

This comparison is shown in Figure 16 between our 
three main types and the same predictions from Benson & 
Devereux (2010). The way to view this comparison is that 
the data has the same line type as in Figure 16, and the pre- 
dictions of the same quantities are shown as double lines. 
The green doubled dashed line shows the evolution of the 
disk fraction, as defined by B/T < 0.5. The red short dashed 
line shows predictions for the evolution of the elliptical frac- 
tion, as defined by B/T > 0.9. The merger fraction predicted 
in the simulation are galaxies defined as having had a major 
merger in the past 1 Gyr, similar to what would be predicted 
for producing a peculiar galaxy (Conselice 2006b). 

The comparison between these quantities is fairly good 
using these definitions, which are somewhat arbitrary as 
we do not know, nor have an obvious way, to compare 
which B/T ratios correlate with early/late type disks. This 
is demonstrated by the white dot-dashed, which is the frac- 
tion of systems with 0.5 < B/T < 0.9, which ranges from a 
fraction of roughly 0.4 to 0.6 within our range of interest. If 
we consider these systems are 'early-type', or based on the 

^ The exact model we use is version 0.9.0.r249 of Galacticus. 
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Figure 16. The comparison of our WFC3 morphologies, probing 
the rest-frame optical at 1 < 2 < 3 and model predictions from 
Benson & Devereux (2010). The single lines are the data, as shown 
in Figure 5, where the dashed line red line is the E/Compact 
systems, the solid blue line denote the peculiar systems and the 
dashed green line is for those galaxies classified as disks. The 
double dashed lines, with the same colour and line type as the 
data, are the model predicts from Benson & Devereux (2010). 
In the models the E/Compact is for systems with B/T > 0.9, 
and the disks are systems with B/T < 0.5, and the peculiars are 
systems which have undergone a merger in the past Gyr. The 
black dot-dashed line shows the fraction of systems with 0.5 < 
B/T < 0.9. 



results in §4.3 as indeed true disks which appear as early- 
types in our classification due to the effects of redshift, then 
the predicted fraction would agree quite well with the early- 
type fraction we observe. However, previous work on more 
detailed comparisons between mergers and theory predic- 
tions show a worse agreement (Bertone & Conselice 2009). 



5 SUMMARY 

This paper presents an analysis of galaxy morphology and 
structure at 1 < z < 3 using WFC3 data utilising several 
popular approaches, include visual estimates of morphology, 
CAS parameters, parametric Sersic profile fitting, as well as 
spectral-type fitting to spectral energy distributions. We fur- 
thermore analyse morphological fe-corrections between ACS 
and WFC3, essentially probing the difference in rest-frame 
morphology between the ultra-violet and optical light at 
z > 1. Our major findings include: 

I. In terms of morphological fc-corrections, a large fraction 
of systems classified as 'peculiar' at 1 < z < 3 within 
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ACS imaging appear as early-types in the WFC3 imaging. 
This result holds to some degree even after convolving the 
ACS image with the larger WFC3 PSF, and by simulating 
nearby galaxies to high-z, and redoing classifications. How- 
ever, these 'early-types' are clearly not the same as those at 
« = in terms of their stellar populations, showing evidence 
for star formation and ongoing assembly. 

II. We further find that these 'early-type'classifications dom- 
inate the non-peculiar galaxy population at 1 < 2 < 3, with 
a fraction near 30%, with a similar fraction of peculiars. The 
remaining systems are classified as 'compact' galaxies, and 
disk-like galaxies, which are both on average 10-20% of the 
population. 

III. We show that there is a significant difference in the 
population of galaxies classified as early /mid/late types as 
defined by spectral-types and those defined through visual 
estimates. The most common difference is found for galaxies 
classified as elliptical - a significant fraction of these have a 
disk spectral type, or even a starburst spectral type. This 
shows that very few purely passive massive galaxy exist at 
1 < 2 < 3, at least within the UDF. 

IV. Examining the CAS parameters, we find that the change 
in CAS values between the 2350 and Hieo bands is similar, 
per unit wavelength, as found between the rest-frame UV 
and optical at 2 < 1 (Conselice et al. 2008). 

V. By using a strict definition of elliptical and disk galaxies 
that matches the concentration, asymmetry, and spectral 
types for nearby systems, we find very few 'formed' systems 
at 2 > 1. We find 10% of our sample would be classified as 
nearby disks at 2 > 1, compared with 75% in the modern 
universe. Likewise, we find < 5% of systems classifiable as 
modern ellipticals, about a factor of three lower than today. 

VI. We find that the inferred merger fraction for the most 
massive galaxies (> 10^" Mq ) at 2 > 1 is between 20-30%. 
This is very similar to previous findings using UDF ACS and 
NICMOS HDF imaging, revealing that even when probing 
the rest-frame optical at 2 > 1 we find a similar merger 
history, and therefore a similar contribution of mergers to- 
wards forming galaxies, as when looking at lower resolution 
and bluer wavelengths. 

VII. We find that the properties of visually classified ellip- 
tical peculiars, and the peculiar ellipticals are very similar, 
suggesting that these galaxies are in some common forma- 
tion scenario. It is very likely that this common process is 
galaxy merging, with each type in a different phase of merg- 
ing. We are therefore likely seeing the formation modes of 
elliptical galaxies through these systems. 

We do not find a large fraction of visually classified 
disk galaxies at 2 > 1 (10-15%). This is potentially due to 
misidentifying ellipticals as rotating disks, however the spiral 
structure, and star forming regions seen in disks at 2 < 1 
and at times at 2 > 1 should be visible within our imaging. 
The number of edge-on disk galaxies, and the high number 
of inclined systems we see, suggests that perhaps we are 
misidentifying some face-on disks as ellipticals, however, the 
mode of formation of these systems must be quite different 
from spirals in the more nearby universe. 



Overall, our conclusion is that visual estimates of mor- 
phology are limited, both because of the PSF of WFC3 limits 
our ability to resolve distant galaxies, as well as the intrinsic 
nature of morphology at higher redshifts. We demonstrate 
that the best way to find distant galax;y 'type' analogs to 
modern galaxies, and therefore to address the question of 
how and when the Hubble sequence forms is through quan- 
titative measurements of both the stellar populations and 
quantitative structures, and leaving behind visual estimates. 
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